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£SJ ■ The two-flavor color superconductivity is examined over a wide range of baryon density with 

a single model. To study the structural change of Cooper pairs, quark correlation in the color 
f*"^ ' superconductor is calculated both in the momentum space and in the coordinate space. At extremely 

high baryon density (~ O(10 10 po)), our model becomes equivalent to the usual perturbative QCD 
treatment and the gap is shown to have a sharp peak near the Fermi surface due to the weak- 
coupling nature of QCD. On the other hand, the gap is a smooth function of the momentum at 
lower densities (~ O(10po)) due to strong color magnetic and electric interactions. The size of the 
Cooper pair is shown to become comparable to the averaged inter-quark distance at low densities, 
which indicates a crossover from BCS to BEC (Bose-Einstein condensation) of tightly bound Cooper 
pairs may take place at low density. 
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Because of the asymptotic freedom and the Debye screening in QCD, deconfmed quark matter is expected to be 
realized for baryon densities much larger than the normal nuclear matter density [jj]. This weak-coupling fermionic 
matter is, however, qualitatively different from a perturbed free Fermi gas system, if the temperature is low enough. 
This is because any attractive quark-quark interaction in the cold quark matter causes an instability of the Fermi 
surface due to the formation of Cooper pairs and leads to the color superconducting phase ||, ||, ^] . 

Current understanding on the color superconductivity has been based on two different theoretical approaches. 
One is an analysis by the Schwinger-Dyson equation with perturbative one gluon exchange, which is valid in weak 
coupling limit at asymptotically high densities Q . The dominant contribution to the formation of Cooper pairs comes 
from the collinear scattering through the long range magnetic gluon exchange q|. In such a weak-coupling regime, 
formation of Cooper pairs takes place only in a small region near the Fermi surface. The other approach is the 
mean-field approximation with the QCD inspired 4- Fermi model introduced to study lower density regionsp, M. In 
this approach, magnitude of the gap becomes as large as 100 MeV and is almost constant in the vicinity of Fermi 
surface. 

The main purpose of this talk is to discuss the superconducting gap over a wide range of baryon density with a single 
model and to make a bridge between high and low density regimes || . To this aim, we make an extensive analysis on 
the structural change in spatial-momentum dependence of a Cooper pair in 2-flavor color superconductivity at zero 
temperature. The momentum dependence of the gap, diffuseness of the Fermi surface, quark-quark correlations in the 
superconductor are the characteristic quantities reflecting the departure from the weak-coupling picture. In particular, 
we show that the spatial size of Cooper pairs which is calculated from the quark-quark correlations, indicates a clear 
deviation from the weak-coupling BCS theory. More complete and detailed analysis was done in Ref. M. 



I. INTRODUCTION 
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II. GAP EQUATION WITH SPATIAL MOMENTUM DEPENDENCE 

Using the standard Nambu-Gor'kov formalism with a two component Dirac spinor \l/ = (tp, %j) ), the quark self-energy 
S(fcp) with the Minkowski 4-momentum k p satisfies 



S (*V>) = ~* / 7^4 9 2 KS{q p )YlDf v (q p - k p ), (1) 



j^k 9 2 KS{ qp )nDf v 

where D?~, is the gluon propagator in medium, S(q p ) is the full quark propagator, and T p is the quark-gluon vertex. 
We study the gap function in the flavor anti-symmetric, color anti-symmetric and J = + channel (the most attractive 
channel within the one-gluon exchange model) j^] : 

A(fc p ) = (A 2 r 2 C 75 ) (A + (fc p )A + (fc) + A_ (k p ) A_ (*)) , (2) 

where t 2 is the Pauli matrix acting on the flavor space, A 2 is a color anti-symmetric Gell-Mann matrix, and C is the 
charge conjugation. A±(k) = (1 ± k • ot)/2 is the projector on positive (+) and negative (— ) energy quarks. 

For the vertex in Eq. ([t|), we use T° = diag.(7 M A a /2, — (7 M A a /2)*). For g 2 in Eq. ([j]), we use a momentum dependent 
coupling g 2 (q, k) in the "improved ladder approximation" || (/3 = (HN C — 2Nf)/'S — 29/3): 
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92 ^ k)=1 ^t\n(( P ^+ P l)/K 2 y ^a X = max(g,fc), (3) 

where p 2 plays a role of a phcnomcnological infrared regulator which prevents the coupling constant from being too 
large at low momentum q,k ~ Aqcd. At high momentum, g 2 shows the same logarithmic behavior as the usual 
running coupling with A identified with Aqcd- The quark propagator in the improved ladder approximation in the 
vacuum is known to have a high momentum behavior consistent with that expected from the renormalization group 
and the operator product expansion. We adopt A=400 MeV and p 2 =1.5 A 2 which are determined to reproduce the 
low energy meson properties for Nf = 2 vacuum fio| . 

The gluon propagator in Eq. (|l|) in the Landau gauge reads, 

pT ph 

D„Jk p ■ k < \k\) = 5 ^ 9 ^ - , (4) 

' 17 k 2 + iM 2 \k \/\k\ k 2 +ml 

where are the transverse and longitudinal projectors. The longitudinal part of the propagator (the electric part) 
has a static screening by the Debye mass m 2 ^ = (Nf /2ir 2 )g 2 fi 2 , while the transverse part (the magnetic part) has 
dynamical screening by the Landau damping M 2 — (7r/4)m 2 ) . This form is a quasi-static approximation of the full 
gluon propagator in the sense that only the leading frequency dependence is considered ITT). 



A. Full gap equation 

The gap equation under the approximations shown above is obtained from the 2-1 element of the Schwinger-Dyson 
(matrix) equation ([!]) which reads 

A±(fc)= [°° dg V±(g, fc;e+ef)^ = A+(<?) 

A, 9 k WE+(q) 2 + \A + (q)\ 2 

r°° A (a) 

Here we used a simplified notation: A±(e^, k) — * A±(k) with E±(q) = q =p /i and being the quasi-particle energy 
as a solution of (e^) 2 = E^q) + A\(e^,q). The explicit form of V± is given in Ref. ||. 

At extremely high density, the Cooper pairing is expected to take place only near the Fermi surface due to weak 
coupling property. In this case, we can safely neglect the antiquark-pole contribution for calculating A + . Furthermore, 
one may replace the momentum dependent coupling constant (||) by that on the Fermi surface g 2 (/i, /i). On the other 
hand, at low densities, sizable diffusion of the Fermi surface occurs and the weak-coupling approximation is not 
justified, and we need to solve the coupled gap equations Eq. (||) numerically. 
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B. Occupation number, correlation function and coherence length 

To clarify the structural change of the color superconductor from high to low densities, it is useful to examine the 
following physical quantities. 

[I] The quark and antiquark occupation number which is related to the diagonal (1-1) element of the quark propagator, 

(t,y)ipf(t,x)) anp er, in the Nambu-Gor'kov formalism; 

"" ( " - K 1 " vm^kw) ■ "+ ,5) = 9( " - *• n - (q) = (6) 

where the superscripts (1,2 and 3) stand for color indices. Since the third axis in the color space is chosen to break 
the color symmetry, quarks with the third color remains ungapped. 

[II] The q-q and q-q correlation functions in momentum space <p±{q) and in coordinate space ip±(r): They reflect the 
internal structure of the Cooper pairs in color superconductor. These correlations are defined through the off-diagonal 
(1-2) element of the quark propagator, (ipf (t, x)ij}j(t, j/)) B uper> 

A ±(<?) .. ,_ s AT f d 3 q A f _ s iqr 



<P±(l) = , <P±(r) = N / — *- 0±(q) e^ r , (7) 

2 v /£±(<z) 2 + |A ± (q)| 2 ' J (2^) 3 * ±W 

where N is a normalization constant determined by J d 3 r\ip+(r)\ 2 = 1. 

[Ill] The coherence length £ c characterizing the typical size of a Cooper pair. It is defined simply as a root mean 
square radius of <p+ (r) : 

= fd*r r 2 |^ + (r)| 2 £dk fc 2 \d<p + (k)/dk\ 2 
^ Jd^ + (rW J Q °°dkk2\v + (k)\ 2 ' U 

It can be shown Q that the quark correlation (p+(r) in the weak-coupling limit behaves as 

p+HooJac^-e^, (9) 

where the Pippard length is given by £ p = (7rA_|_(/x)) . 

In a typical type-I superconductor in metals, the Pippard length is of the semi-macroscopic order £ p ~ 10^ 4 cm, 
whereas inverse Fermi momentum is of the microscopic order k^ 1 ~ 10 -8 cm. Besides, there is another scale wq, the 
Debye cutoff, which limits the range of attractive interactions to be just around the Fermi surface |A; — &f| < Wd- 
The inverse of the Debye cutoff is in between the two scales above: Wp 1 ~ 10~ 6 cm. Therefore there is a clear 
scale hierarchy, A -C lo-q <C fcp. On the other hand, since there is no intrinsic scale lo-q in QCD, scale hierarchy 
at extremely high density simply reads A ~ [ie~ c l 9 <C kp ~ [i. At lower densities, however, such scale separation 
becomes questionable for g is not small. 



III. MOMENTUM DEPENDENT GAP FROM HIGH TO LOW DENSITIES 

In Fig. |l|(a), we show the gap A + (fe) as a solution of the full gap equation for a wide range of densities. Since the 
actual position of the Fermi surface moves as we vary the density, we use kj \i as a horizontal axis in the figure which 
helps us to understand the change of global behavior. The figure shows that the sharp peak at high density gradually 
gets broadened and simultaneously the magnitude of the gap increases as we decrease the density. 

The characteristic features at high density are (i) there is a sharp peak at the Fermi surface, and (ii) the gap decays 
rapidly but is nonzero for momentum far away from the Fermi surface. The property (i) is similar to the standard BCS 
superconductivity but (ii) is not, due to the absence of intrinsic ultraviolet Debye-cutoff of the gluonic interaction in 
QCD. As for the magnitude of the gap at high density, if one estimates the several contributions to the kernel V± in 
the gap equation (^) separately, one finds that the color-electric interaction enhances the gap considerably. This may 
be understood as follows: In the coordinate space, the Debye-screened electric interaction is a Yukawa potential. Such 
a short-range interaction alone can form only a loosely bound Cooper pairs and a very small gap. However, if the 
magnetic and electric interactions coexist, small size Cooper pairs are formed primarily by the long-range magnetic 
interaction. Then, even the short-range electric interaction becomes effective to generate further attraction between 
the quarks. 
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FIG. 1: (a) A+(fc) as a function of k/u for various densities (i = 2™ A with n = 1, 2, 3, 12. (A =400MeV). All the calculations 
are done with the momentum dependent vertex g 2 (q,k) and with the antiquark-pole contribution, (b) Chemical potential 
dependence of the gap A+(fc = jj.) in the full calculation compared with the analytic result which is normalized at the highest 
density fj, = 2 12 A. 



The characteristic features at low density are that (iii) the sharp peak at the Fermi surface disappears, and (iv) all 
the contributions neglected in the weak-coupling limit are actually not negligible. A close look at the effects of each 
contributions to the gap leads us to a conclusion that the color superconductivity at low density is not a phenomenon 
just around the Fermi surface. This is confirmed by computing the occupation number following eq. (^|). The result 
shows that the Fermi surface is diffuse substantially at low density. 

In Fig. |l|(b), the gap at the Fermi surface A + (/i) is shown as a function of the chemical potential. It decreases 
monotonically as \i increases as is shown on the figure, but turns into an increase at much higher density /i > 10 6 
MeV. An analytic solution based on the weak coupling approximation (includes only the magnetic gluons, fixed 
coupling, and ignores the anti-quark pole contribution) is also shown in Fig. Mb). The magnitude of the analytic 
solution is normalized to the numerical solution at the highest density /i = 2 A ~ 1.6 x 10 6 MeV. At high density, 
/i-dependence of the numerical result is in good agreement with the analytic form which has a parametric dependence 
A + (/i) oc <?~ 5 /i exp(— 3n 2 /\/2g) with g 2 = g 2 {p, fS). On the other hand, the difference of the two curves at low density 
implies the failure of the weak-coupling approximation. 



IV. CORRELATION FUNCTION AND COHERENCE LENGTH 



One of the advantages of treating the momentum dependent gap is that we are able to calculate the correlation 
function which physically corresponds to the "wavefunction" of the Cooper pair. Such correlations have been first 
studied in Ref . Jl2| ] in the context of the color superconductivity (but with much simpler model than eq. (^J) , and only 
at lower density) - ^ 

Figure §(a) shows the coherence length £ c of a quark Cooper pair denned as the root mean square radius of the 
correlation function [see Eq. (0)] . The size of a Cooper pair becomes smaller as we go to lower densities. This tendency 
is understood by the behavior of the Pippard length £ p — lj_n A + (/i) (which gives a rough estimate of the coherence 
length) together with the behavior of A + (/i) shown in Fig. Mb). 

Also, the Cooper pair becomes smaller as density increases beyond /i = 2 12 A. However, it does not necessarily 
imply the existence of tightly bound Cooper pairs. In fact, the size of a Cooper pair makes sense only in comparison 
to the typical length scale of the system, namely the averaged inter-quark distance d q for free quarks denned as 



a \ 1/3 



As we go to higher densities, the ratio S, c /d q increases monotonically as shown in Fig. ||(b). Namely loosely bound 
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FIG. 2: (a): Density dependence of the coherence length, (b): Ratio of the coherence length £ c and the average inter quark 
distance d q as a function of the chemical potential. 



Cooper pairs similar to the BCS superconductivity in metals are formed at extremely high densities. (Recall that the 
typical ratio in superconductivity is fcp/A ~ 10 4 ) 

At the lowest density in Fig. 0, the size of the Cooper pair is less than 4 fm and the ratio £, c /d q is less than 10. This 
is not similar to the usual BCS system. The transition from (, c /d q 3> 1 to £ c /gL ~ 1 as M decreases is analogous to 
the crossover from the BCS-type superconductor to the Bose-Einstein condensation (BEC) of tightly bound Cooper 
pairs |n| . Our result here suggests that the quark matter possibly realized in the core of neutron stars may be rather 
like the BEC of tightly bound Cooper pairs. 

For better understanding of the internal structure of the quark Cooper pair, let us consider the correlation function 
in the coordinate space. Fig. [| shows the spatial correlation of a Cooper pair at various chemical potentials normalized 
as J d 3 r\ip + (r)\ 2 = 1. 
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FIG. 3: The quark-quark correlation function <^+( r ) m the coordinate space (p+(r) for several different chemical potentials. 
<p+(r) is normalized to be unity. 



At high density most of the quarks participating in forming a Cooper pair have the Fermi momentum fcp = A 4 
giving a sharp peak in the momentum space correlation. In the coordinate space, this corresponds to an oscillatory 
distribution with a wavelength A = 1/fj, without much structure near the origin. (The oscillation is also evident 



6 



from the factor sin(/nr) in the approximate correlation function Eq. (H).) At lower densities, accumulation of the 
correlation near the origin in the coordinate space is much more prominent in Fig. ||. This implies a localized Cooper 
pair composed of quarks with various momentum. 



V. SUMMARY AND DISCUSSION 



We have studied the spatial- momentum dependence of a superconducting gap and the structure of the Cooper pairs 
in two-flavor color superconductivity, using a single model for a very wide region of density. Nontrivial momentum 
dependence of the gap manifests itself at low densities, where relatively large QCD coupling allows the Cooper pairing 
to take place in a wide region around the Fermi surface. Note that our results can be easily extended to Nf = 3 
case. Our results imply that the quark matter which might exist in the core of neutron stars or in the quark stars 
could be rather different from that expected from the weak-coupling BCS picture and could be more like a BEC of 
tightly bound Cooper pairs. Study of the finite T phase transition of this strong coupling system is currently under 
investigation. 
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